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ABSTRACT 

The standard model of planet formation considers an initial phase in which planetesimals form from 
a dust disk, followed by a phase of mutual planetesimal-planetesimal collisions, leading eventually 
to the formation of planetary embryos. However, there is a potential transition phase (which we 
call the "snowball phase"), between the formation of the first planetesimals and the onset of mutual 
collisions amongst them, which has often been either ignored or underestimated in previous studies. 
In this snowball phase, isolated planetesimals move on Kcplerian orbits and grow solely via the direct 
accretion of sub-cm sized dust entrained with the gas in the protoplanctary disk. Using a simplified 
model in which planetesimals arc progressively produced from the dust, we consider the expected sizes 
to which the planetesimals can grow before mutual collisions commence and derive the dependence of 
this size on a number of critical parameters, including the degree of disk turbulence, the planetesimal 
size at birth and the rate of planetesimal creation. For systems in which turbulence is weak and the 
planetesimals are created at a low rate and with relatively small birth size, we show that the snowball 
growth phase can be very important, allowing planetesimals to grow by a factor of 10 6 in mass before 
mutual collisions take over. In such cases, the snowball growth phase can be the dominant mode 
to transfer mass from the dust to planetesimals. Moreover, such growth can take place within the 
typical lifetime of a protoplanetary gas disk. A noteworthy result is that, for a wide range of physically 
reasonable parameters, mutual collisions between planetesimals become significant when they reach 
sizes ~ 100 km, irrespective of their birth size. This could provide an alternative explanation for the 
turnover point in the size distribution of the present day asteroid belt. For the specific case of close 
binaries such as a Centauri, the role of snowball growth could be even more important. Indeed, it 
provides a safe way for bodies to grow through the problematic ~1 to 50 km size range for which the 
perturbed environment of the binary can prevent mutual accretion of planetesimals. From a more 
general perspective, these preliminary results suggest that an efficient snowball growth phase provides 
a large amount of "room at the bottom" for theories of planet formation. 
Subject headings: planets and satellites: formation 



1. INTRODUCTION 

The standard core-accretion model general ly treats 
plane t formation as occurring in two s tages (jLissauerl 
[1991 IChambersI l200i lArmitagel I20II : In the first 
stage, mountain-sized planetesimals form from small 
dust grains embedded in a gas-rich protop l anetary 
disk (lWeidenschillindfT997tlBlum fe Wurm| [200l lYoudinl 
120081 iChiang fe Youdinl I2009T ). and then in the sec- 
ond stage these planetesimals go on to accrete one an- 
other to form planetary em bryos (jKokubo fc Idal 119961 . 
Il998t iWeidenschillind I1997D which eventually go on to 
merge into full planets dChambers fc Wetherilll 119981 : 
iLevison fc Agnorl l2003t iKokubo et all 12006ft . The first 
stage - planetesimal formation - is of crucial importance, 
as it sets the initial conditions upon which subsequent 
stages of evolution depend. However, the details of plan- 
etesimal formation are still poorly understood and re- 
main somewhat controversial. 

One model for the formation of planetesimals en- 
visages them growing via mutual sticking collision s 
(|Weidenschilling fc Cuzzil Il99l iWeidenschillingl [l997h . 
Laboratory experiments show that ^m sized dust grains 
can stick together efficiently through various surface 
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forces, including the van der Waals and electrostatic 
forces, causing them t o form mm to cm-sized aggregates 
urn fc Wurm!l2CJ0l . Beyond this mm to cm range, 
particles become less sticky while their gravitational in- 
teractions remain w eak, leading to collisi onal disruption 
rather than growth (jDominik et al.ll2007T ). Furthermore, 
cm to m sized objects begin to decouple from the gas, 
so that they experience collisions with higher relative ve- 
locities and also begin to experience significant aerody- 
namic drag which causes them to quickly spiral into the 
protostar on a timesca le of a few hundred orbital periods 
(| Weidenschillingll 1 9 77| ) - the well-known "meter-barrier" . 

An alternative scenario that has been suggested to 
try and circumvent this meter-barrier, envisages that 
the km-sized planetesimals form directly via gravita- 
tional instabilities in a dense particle s ub-disk near the 
mid-plane of the protoplanetary disk (|Safronovl 119691 : 
IColdreich fc "Wa7^rT978l : lYoudin fc Shull2002D . However, 
it has been pointed out that even low levels of tur- 
bulence would prevent solids fro m settling down to a 
sufficiently dense mid-plane la yer (I Weiden schilling] 1 1 9801 : 
iCuzzi fc Weidenschillind" 120061 ). 

Recently breakthroughs have been m ade in both sce- 
narios . For the collisional scenario, iTeiser fc Wurml 
(|2009f ) found in experiments that high-velocity (tens of 
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m.s -1 ) collisions between small dust particles (< cm) 
and a pre-planetesimal body (> dm) can lead to efficient 
growth, suggesting that some "lucky" pre-planetesimals 
may escape from ero sive collisions and go on to form km- 
sized planetesimals (jJohansen et aljT2008[) . 

For the instability scenario, it is found that turbu- 
lence can act to concentrate dust, thus providing a high 
density, low velocity dispersion environment in which 
planetesimals may form via various mod els of instability 
(jJohansen et al.l 120071 ICuzzi et all 120081 ). The planetes- 
imals formed in these new turbulent models can be 1-2 
orders of magnitude larger than those formed in tradi- 
tional models. 

Planctcsimal formation is thus far from being com- 
pletely understood and current models arc a work-in- 
progress. However, a point worth mentioning is that 
none of these different scenarios ensure that all plan- 
etesimals appear at the same time at a given location. 
In fact, any model which requires non-gl obal density en- 
hance ments to form planetesim als, e.g. Uohansen et al.l 
(|2007t) and ICuzzi et al.l (|2008[) . implicitly assumes that 
certain regions of the disk will form planetesimals first, 
while the rest of the d isk remains "dust y" . F urthermore, 
two re cent studies by IChambersI (j2010D and ICuzzi et"aLl 
(l2010h argue that, for the turbulent concentration sce- 
nario, planet formation efficiency (set by the probabil- 
ity for turbulently-formed clumps of mm-sized grains to 
collapse into planetesimals) could be relatively low, de- 
pending on the values of several crucial parameters such 
as the dust-to-gas density ratio, disk viscosity, density 
and profile. This means that, at a given location in the 
disk, there could be a wide time gap between the moment 
when the first planetesimals form and the moment when 
most of the available solid mass has been converted into 
planetesimals. 

There could thus be a long transition period dur- 
ing which individual planetesimals are embedded in a 
disk where most of the mass of solids is still in small 
grains. This issue has tended to be ignored in stud- 
ies investigating the next stage of planet formation, 
i.e., planetesimal accretion. Most of these studies con- 
sider a system where all planetesimals are present at 
time to, possibly with a distribution of initial sizes, 
which then grow by pairwise accr e tion (iWeidenschillina 
[19971 IKokubo fe Idal [19961 [l998t IBarnes et al.1 12009ft . 
Nevertheless, ther e are several noteworthy exceptions. 
iWetherill fe Inabal (j2000D considered that planetesimals 
progressively appear over a 10 5 years timescale, but only 
considered mutual planetesimal accretion as a possible 
growt h mode, implicitl y negle cting the contribution of 
dust. IMorbidelli et al.l (|2009D . while attempting to fit 
the known asteroidal size-distribution, did consider in 
one of their simulations the possible accretion of dust by 
planetesimals which formed early, but did so only for one 
specific case: large, 100 km size d seed-planetesimals pro- 
duced over 2 Myr. In addition, iLeinhardt fe Richardson! 
(|2005[) attempted to model accretion of dust onto plan- 
etesimals as well as planetesimal-planetesimal collisions 
but they started with very large planetesimals. The 
most promising study has been recently performed by 
IPaardekooper fc Leinha rdt (2010]), who envisage dust ac- 
cretion onto planetesimals more explicitly, but restricted 
to the specific case of a close-in binary, and using only a 



simplified 2-dimensional model. 

In this study, we plan to take these pioneering studies a 
step further and investigate in detail planetesimal growth 
during the transition period when isolated planetesimals 
and a primordial dust-disk coexist. We outline a new 
growth mode, which dominates during this early phase, 
during which planetesimals experience negligible gravi- 
tational or collisional interactions with one another, and 
grow mainly (or solely) via the accretion of dust or ice 
that they sweep up— in the manner of a rolling snowball. 
In this paper, we refer to this dust-fueled planetesimal 
growth phase as the "snowball" growth phase. 

The paper is organized as follows. The snowball 
growth rate (growth only by dust accretion) is derived, 
using a semi-analytical approach, in §2. Next, in §3, 
we show how long the snowball phase could last and to 
what radial size and mass fraction the planetesimals can 
grow via this snowball growth mode. Then, in § 4, the 
implications for planet formation in both the Solar Sys- 
tem and in close binary systems are discussed. Finally, 
we conclude in § 5. Additionally, all variables defined in 
this paper are listed alphabetically in table [1] 

2. MODEL 
2.1. Planetesimal Formation Rate 

We adopt a simplified analytical model which contains 
only two components: dust and planetesimals. The dust 
in our model is assumed to be in the sub-mm to m m size 
range , as (1) this is the typical size of chondrules (jScottl 
l2007h . and (2) particles of this approximate size settle 
down to the mid-plane of the disc on a relative short 
timescale of ~ 10 3 yr, and can survi ve the inward drift 
for as long as - 10 6 yr (|Youdinll2008T) . 

Individual planetesimals are assumed to be quickly 
produced from du st, on a timescale it And, whether 
thro ugh collisions (iTeiser fc Wurml I2009T) or instabili- 
ties (j Johansen et al.l 120071 : ICuzzi et al. I l2008h . A rea- 
sonable referen ce value for tj.ind might be ~ 10 4 yr 
(jLissauerl 119931 ). However, we do not restrict our study 
to tf.ind = 10 4 yr. In fact, as we show in section 3, the 
precise value of tf.ind is unimportant to the final results 
as long as it is not too long (i.e. < t snow , see §3). 

We make the simplifying assumption that, once formed 
from the dust, all planetesimals have the same size R p o, 
which wc take as the initial condition for newly-formed 
planetesimals in our model. To account for the dispersion 
in times at which planetesimals are created in the sys- 
tem, we introduce an efficiency factor e p (0 < e p < 1) for 
planctcsimal formation at a given location in the proto- 
planctary disk. We then consider a given radial distance 
in the disk and define the local planetesimal formation 
rate as 



dN 
~dt 



pO \^f.ind 



M, 



(1) 



where N and M p q are the surface number density and 
initial individual mass of the planetesimals respectively, 
and Erf is the du st surface density. Further, following 
I Chambers] (|2010l) . tf is the characteristic planetesimal 
formation timescale defined as 
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The surface number density of planctesimals in the sys- 
tem thus increases linearly according to the relation 



N 



t 



pO \ 6 f.ind 



pO 



(3) 



assuming that N = at t = 0. 

Note that t / does not represent the time it takes for an 
individual planetesimal to form (that is t/.j nc j) but rather 
the time it takes for most of the dust to be converted 
into planetesimals, as parameterized by the factor e p . If 
€. p = 1, i.e., tf = t f.i n d, this implies that all the plan- 
etesimals are created at the same instant. This should 
be treated as an idealized or limiting case, and in real - 
ity, we expect e p < 1, i.e., tf > tf. in d- iChambersI (|201QO 
finds that tf can easily vary over 10 orders of magnitude, 
from 10 4 to 10 14 yr, depending on local conditions in the 
protoplanetary disk. 

Our model, assuming a constant rate of planetesi- 
mal creation as well as a single initial planetesimal size, 
is very simplified. However, it is accurate enough for 
the order-of-magnitude approach adopted in this present 
study and allows us to identify and quantify the snowball 
phase in a convenient manner. 

2.2. Growth Rate via Dust Sweeping 

When a planetesimal is sweeping through a dust disk, 
the mass growth rate through dust accretion is 



AL = E d -KRi 



2H, 



■Vrel, 



(4) 



where M p , and R p arc the mass and radii of the plan- 
etesimal, Hd is the scale height of the dust-disk, and v re i 
is the relative velocity between planetesimal and dust. 
Ed is a dimensionlcss coefficient accounting for the effi- 
ciency of dust accretio n onto planetesimals. According to 
iTeiser fc WuTrnl (f2009h . accretion can be very efficient as 
long as the dust size is smaller than mm. Hence, through- 
out this paper, we take a medium value Ed = 0.5. From 
Eqn|4] we can also derive the growth rate in the plan- 
etesimal radius as 



R p = 



Ed 



El 



Vrel, 



(5) 



where is the planetesimal density and H g is the scale 
height of the gas-disk. Throughout this paper, we adopt 
/)* = 3g.cm -3 and use a protoplanetary dis k scaled 
by the Minimum Mass Solar Nebula (MMSN, lHavashil 
dU)), which gives S d ~ 10/ d / Jce (a/AU)- 3 / 2 g.cm- 2 , 
and H g ~ 0.05(a/AU) 5/4 AU, where f d is the scaling 
factor of the disk mass relative to the MMSN, and fi ce 
accounts for the enhancement of solid density beyond 
the ice line, ai ce . We take L nP . = 4.2 if a > ai ce , and 
f lce = 1 otherwise (|Ida fc Linll2004| ). Note that in EqnH 
and Eqn(5l we do not consider the gravitational focusing 
effect, since dust is well coupled with the gas while the 
focusing effect is only well applied to 2-body dynamics 
in which only the gravity of the two objects themselves 
matters. As a consequence, v re i is basically the local 
differential velocity between a large planetesimal and a 
gas streamline. Assuming an axisymmetric and pressure- 
supported gas disk, this leads to 



V re l 



(H g /a) 2 v k ~ 75m.s" 



(G) 



where Vk is the local Keplcrian velocity at a AU. Substi- 
tutinging this into EqnJSl we get the growth rate in the 
planetesimals' radii due to the direct accretion of dust, 



R P ~7xW- 7 fdf lce (J^j 



-11/4 



km.yr . (7) 



Taking f d x f lce - 1 - 10 and H g /H d ~ 1 - 1000, we get 
Rp ~ 10~ 6 - lO^km.yr" 1 at 1 AU. 

It is worth noting from EqnJT] that the planetesimal 
radii follow a linear growth rate. Since the planetesimal 
formation rate is also assumed linear as shown in Eqn(31 
then we can derive the mass-weighted average planetesi- 
mal radius as (see Appendix for detail) 



(Rp) ~ R 



p0 



1/3 



Rvt, 



(8) 



as well as the Planetesimal Mass Fraction (PMF) with 
respect to the total solid mass, 
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(Rp) and PMF are two key statistics which trace the 
typical size and total mass of planetesimals in the proto- 
planetary disk. Note that, for the growth rate we derive 
in EqnJTl a constant dust surface density is implicitly as- 
sumed. In reality, dust would be consumed both through 
the production of new planetesimals and through the 
growth of planetesimals which formed earlier. To account 
for this depletion, we use PMF as an alarm or flag: if 
PMF > 0.5, we then understand that a large fraction 
of dust has been consumed and we turn off planetesimal 
growth via dust accretion. 

2.3. Growth via Mutual (Planetesimal-Planetesimal) 
Accretion 

As planetesimals begin to populate the disk, they can 
begin having mutual encounters. The timescale for the 
onset of mutual collisions decreases as the number den- 
sity increases. Assuming equipartition between in-plane 
and out-of-plane motions, we get 



i coZ ~l/ {nAVTr(Rp) 2 ) 

2 / Jrl'^nP* \ ( t 

~ 3tt 
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(—)' 
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yr, (10) 



where n ~ N/(2ai p ) and AV ~ 2i p Vk are the volume 
number density and average relative velocity of the plan- 
etesimals, respectively, and i p is their average orbital in- 
clination. 

3. SNOWBALL GROWTH 

We define snowball growth as beginning with the for- 
mation of the first planetesimal and ending at the point 
when mutual planetesimal collisions take over as the 
dominant growth mode, i.e., when t = t co i. If the snow- 
ball growth is efficient, i.e., (Rp) S> R P o, then the dura- 
tion of the snowball phase t s 
derived as 



can be approximately 
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Given t snow and using Eqn|S] and EqnJHl the mass- 
weighted size (Rsnow) and the mass fraction of planetesi- 
mals (PMF snow ) at the end of the snowball growth phase 
can be estimated as, 

Rsnow ~R P o + 0.7 (/ d f lce y 1/4 (j^) (ten) 

j a r y /2 m 3/4 kmi (12) 

l lO-'Wyr- 1 I V AU/ ' v ' 

and 

PMf._~2.4x 10- (&) (JL) V4 ,13) 

We wish to emphasise that: 

1. EqnJTTJ EanfT2l and Ecmll3l are approximate solu- 
tions which are valid only if (R p ) ^> R p0 . When 
this approximation is invalid ((Rp) ~ R P o), then 
tsnow, Rsnow and PMF 3now cannot be expressed 
analytically and numerical solutions are required 
(see FigHin §3.3). 

2. EqnfTT] and Eqn[T2] depend on tf (not tf.ind). 
The coefficient of 10 4 yr appearing in Eqn|TI] and 
EqnfT2]is not related to the value tf.ind ~ 10 4 yr 
mentioned in section 2.1. In fact, tf.ind does not 
affect the results (Eqn|TIJ EanlTSIand EqnfT3|) pro- 
vided that tf.ind < tsnow- 

3. PMF snow is the total PMF at the end of snow- 
ball phase, i.e. t = t snow , it is not just the PMF 
contributed purely by the snowball growth phase. 
(See also in § 3.1 and § 4.1.1 for details) 

4. Since R p is proportional to /<j and fi ce , R sn ow ac- 
tually increases with the disk density. 

3.1. Efficiency of Snowball Growth 

An efficient snowball growth phase would imply that 
a significant proportion of the disk's solid-mass is con- 
verted into planetesimals in this phase, so one essen- 
tial condition is that PMF snow should be close to unity. 
Let's define this high.-PMF sn0 w condition as correspond- 
ing to 0.1 < PMF S now < 1- From EqfJ31 we see that 
high PMF S now fractions are favored by high values of 
Hg/Hd- At a = 1 AU from the star, for instance, 
0.1 < PMFsnow requires H g /Hd > 40. To a first or- 
der, Hg/Hd is an indicator of the degree of turbulence in 
the protoplanetary disk; the larger H g /Hd, the weaker 
the turbulence. Therefore, snowball growth is more effi- 
cient for a disk that has weaker turbulence (higher value 
of Hg/Hd). This can be simply understood as being the 
condition that, for the same amount of dust mass, this 
mass is concentrated in a thinner, denser disk and is thus 
more easily accreted by planetesimals that always stay 
very close to the mid-plane. 

However, while the above condition is necessary, it is 
not sufficient. This is because, according to Eqn[9] (re- 
placing R p by R snow and t by t snow ), PMF snow has 



contributions from two sources: the snowball growth 
term (R sn ow / Rpo) 3 , and the initial planctcsimal creation 
term t snow /tf. Therefore, the efficiency of snowball 
growth should be measured using both PMF sn0 w and 
the ratio R sn ow/RpO- In this paper, we define the crite- 
ria for an efficient snowball growth phase as (1) 0.1 < 
PMFsnow < 1-0 (or equivalently, H g /H d > 40), and 
(2) Rsnow/RpO > 10 (or equivalently, t snow /t f < 10~ 3 ). 
Note that although the efficiency of snowball growth 
decreases with the ratio of t snow /tf, the efficiency ac- 
tually increases with the absolute value of t snow , since 

Rsnow OC tsnow 

3.2. Examples 

As shown in Eqn ll3l and discussed in the above sub- 
section, the efficiency of the snowball growth phase de- 
pends on the ratio between the gas and dust scale heights 
(Hg/Hd) which is an indicator of the degree of disk 
turbulence. Here we consider three cases with weak 
(Hg/Hd — 150), intermediate (H g /Hd = 15) and strong 
(Hg/Hd = 1-5) disk turbulence, which lead respectively 

to Rp = 10" 4 , 10" 5 and 10~ 6 km.yr" 1 at 1 AU in a 1 
MMSN disk. For each case, we then consider two further 
sub-cases with different characteristic planetesimal for- 
mation timescales (not tf.ind)- tf = 10 5 yr and tf = 10 11 
yr and with the same initial planetesimal size of R p q = 1 
km. Solving for t co i, (Rp) and PMF as an explicit func- 
tion of time (using EqnJHl Eqn[9] and Eqn ITO]) . we plot 
the results in Fig[TJ 

In Figdl the locations of the triangles (tf = 10 5 yr) and 
squares (tf = 10 11 yr) mark the important point at which 
the system transits from snowball growth to the mutual 
collision mode. As expected, the snowball phase is much 
more efficient in the weakly turbulent (high H g /Hd) case. 
Additionally, if the initial planetesimal creation process 
is inefficient (tf = 10 11 yr), then the snowball phase is 
longer and planetesimals can grow up to R sn ow ~ 40 
km purely via dust accretion, whereas when the creation 
process is efficient and tf = 10 5 yr, the planetesimals can 
only reach ~ 3 km before mutual collisions take over. At 
the end of the snowball growth phase, both cases have 
comparable values of PMF snow (~ 35% and ~ 45% re- 
spectively), both of which are below our alarm value of 
PMF = 0.5. 

Conversely, for the highly turbulent case, snowball 
growth is less efficient. Although the snowball phase 
lasts longer than in the low turbulence case, planetes- 
imals don't grow quickly enough to reach sizes bigger 
than ~ 1.2 km (tf = 10 5 yr case) or ~ 5 km (tf = 10 11 
yr case). 

In summary, Fig[T]demonstrates the general trend that 
snowball growth is more efficient when (a) planetesimal 
creation rates are lower (i.e., larger tf), and (b) when 
disk turbulence is weaker (i.e., larger H g /Hd). In addi- 
tion, FigjT] provides a method (by finding the crossing 
point in the t - t col plane) by which t snow , R sn ow, and 
PM Fgnow can be found more accurately than just us- 
ing the approximate expressions of Eqn llU Eqn ll2l and 
Eqn[ni 

3.3. Mapping t snow and R sn ow 

The results of § 13.21 were obtained for a fixed initial 
planetesimal size of 1km. We now extend these results 
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to a more general case in which the initial planetesimal 
size, Rpo, is allowed to vary. Such a consideration is im- 
portant for two reasons: Firstly because R sn0 w does not 
vary linearly with R p o-, and secondly because the initial 
planetesimal size is a poorly constrained parameter that 
strongly varies from one planetesimal formation scenario 
to the other. 

For illustrative purposes, we focus here on the most 
favorable case for snowball growth, i.e., weak turbulence 
{Hg/ Hd = 150), and vary both R p q and tf as free param- 
eters. FigJ3] shows t snow and R sn0 w in the R p q — tf plane 
and was generated by numerically solving the equation 
t = t co i. As the values of t snow and R sn0 w plotted in 
Figj2] are accurately solved using this numerical method, 
they do not rely on the approximation (R p ) 3> R p o 
which was assumed in Eqn[TT] and EqiifT^l Note that 
in the bottom-right corner of Fig(2J the R sn ow contours 
become vertical, while the t snow contours become hori- 
zontal. This occurs because we suppress snowball growth 
whenever the PMF becomes greater than the alarm value 

O.t).t sr iow ^ ttran ^° 10 10 

Not surprisingly, snowball growth is at its most effi- 
cient in the top part of the figure, where tf is long (or 
cquivalently, the planetesimal creation rate is low). In 
addition, from EqnfT2]we see that R S now/R P o °c i?~ ' 25 , 
so that snowball growth, measured in terms of the in- 
crease in planetesimal size, is more efficient for smaller 
initial planetesimals, i.e. on the left-hand side of the 
plot. This result should be expected since, for a given 
total mass of solids, smaller initial sizes lead to larger 
total collisional cross-section for dust accretion. 

In contrast, snowball growth gets less efficient when 
starting from larger initial planetesimals. However, even 
for R p q = 100 km, a factor 10 increase in size (1000 in 
mass) is possible for low-turbulence disks with a low 
planetesimal formation rate (tf > 10 11 yr). Setting 
this factor 10 size increase as the criteria for an efficient 
snowball phase (see Sec J3.1[) . we see that the correspond- 
ing minimum value for tf increases from ~ 10 yr for 
R p o = 0.1km to - 10 13 yr for R pQ = 1000 km. 

However, another crucial constraint is that t snow 
should not exceed the time for disk dispersal (probably 
< 10 7 yr). This constraint is relevant for the following 
considerations: if gas is dissipated, then (1) both dust 
and planetesimals would have Keplerian velocities, and 
there would thus be far less efficient dust-sweeping and 
no snowball growth for planetesimals, and (2) one would 
like to be able to form gas giants such as Jupiter in the 
Solar System, and hence require that massive solid plane- 
tary cores be able to form before disk dissipation. Given 
this additional constraint, all the solutions in the top- 
right part of FigJ5]are ruled out. A factor 10 size increase 
in less than 10 7 yr is thus only possible for R p o < 100 km. 

4. DISCUSSION 

4.1. Implications for Planet Formation in Single-Star 
Systems: The Solar System and Extrasolar 
Systems 

We emphasize that this current section of discussion 
(SHI) is based on the results of FigJ5] for the weak tur- 
bulence case where snowball growth is most favorable. 

4.1.1. Dust to Planetesimal Transition Timescale, ttran 



It is worth noting that in FigJ5] we always have 
PMF snow ~ 35% — 50%. As a consequence, in this case 
tsnow can be interpreted as corresponding approximately 
to the typical timescale, ttran, for which the solid disk 
goes from being dust-dominated to being planetesimal- 
dominated. This timescale is especially relevant for 
planet formation scenarios because it can be indepen- 
dently mea sured, through bo th protoplanetary disk ob- 
servations dNatta et al.ll2007t) and meteoritic evidence in 
our S olar System rtCuzzi fc Weidenschilhngl [20061 IScottl 
20075). All these observational studies seem to agree on 
a transition timescale somewhere between 1 and lOMyr, 
i.e., t snow <~ ttran ~ 10 6 — 10 7 yr, which is shown as the 
red region in Figj2] 

Further, we note that if i? srlo , 



R p o then 



ttran ^ t snow ~ tf , 



whereas if R sn0 w 3> Rpo then 



ts 



These relations occur because snowball growth implicitly 
includes two simultaneous processes, i.e. the formation of 
new planetesimals, plus the snowball growth of new plan- 
etesimals. As discussed in §3.1, snowball growth domi- 
nates the snowball phase only if R snow ^ RpQ. This 
requires R p o < 100km for the red region of FigEl 
In contrast to the above, if PMF snou , <C 0.1 then 



tf ^> t snc 



Taking into account all these considerations, one should 
take care not to confuse tf with either t snow or t tra n, as 
it may potentially be orders of magnitude different from 
either. 

Irrespective of the precise scenario in question, we note 
that snowball growth provides an alternative channel 
through which dust can be converted into planetesimals. 
Moreover this can be an efficient process; even if the ini- 
tial planetesimal creation is very inefficient, one can still 
convert most of the solid mass into planetesimal within 
1-10 Myr. 

4.1.2. Planetesimal Birth-Size 

Following the previous chapter, we take as a reference 
the observationally derived constraint of t tra n ~ 10 6 — 10 7 
and the additional result that, for the weak turbulence 
case, t snow ~ ttran- As a consequence, the most prob- 
able location for the dust-to-planetcsimal transition is 
the red region of Fig|2] An important feature is that 
almost the whole of this red region corresponds to a 
100 < Rsnoiu < 1000 km range. This implies that, for the 
weak turbulent case, mutual planetesimal-planetesimal 
collisions (or planetesimal accretion) start when plan- 
etesimals are 100-1000 km in radii, regardless of their 
initial size R p q. 

Interestingly, this impli cation is consistent with the re- 
sult of a recent study by iMorbidelli et al.l (|2009() , which 
found that the size-distribution of asteroids in the So- 
lar System, in particular the knee in this distribution 
around 100 km, could be well reproduced if planetesi- 
mal accretion (and mutual fragmentation) starts from 
100-1000 km-siz e d obj ects. This was interpreted by 
IMorbidelli et~aTl (l2009h as being an indication that as- 
teroids were born big, i.e., they emerge from their dust- 
to-pebbles progenitors with a size exceeding 100 km. We 
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argue that our results could provide an alternative expla- 
nation, i.e., that 100 to 1000 km is the characteristic size 
for the end of snowball growth and the onset of mutual 
collisions, a size that might be much larger than the size 
R p o of the initial seed planetesimals. In other words, the 
results of Morbidclli ct al. (2009) could actually provide 
some supporting evidence for the snowball growth, since 
the result R sn ow = 100 — 1000 km is robust when apply- 
ing the snowball concept to parameters consistent with 
the Solar System. However, this inference depends on 
whether the weak turbulence case is consistent with the 
Solar System. 

As mentioned in §1, a form of sno wball growth was 
taken into account in the simulations of lMorbidelli et all 
(|2009t) (see their Fig. 6b). However, there are impor- 
tant differences between their work and this present 
study. First of all, dust sweeping was only considered for 
large (> 100 km) seed planetesimals, implicitly neglect- 
ing its possible effect on smaller initial bodies. Another, 
perhaps more important difference lie s in the respec- 
tive m odels adopted for dust sweeping. iMorbidelli et al.l 
(l200l argue that it proceeds in a fast "runaway" mode 
- thus making it end much earlier - because the rel- 
ative velocity between the planetesimals and the dust 
(vvel ~ 75ms -1 ) is smaller than the planetesimals' es- 
cape velocity (v esc ~ 100ms -1 for a 100 km size), and 
hence the dust is gravitationally focused onto the larger 
planetesimals. In this study we ignore runaway growth, 
even for large planetesimals, because we implicitly as- 
sume that the coupling of the dust to the gas is strong 
enough to prevent solid grains from being gravitationally 
deflected onto the planetesimals according to the usual 
gas-free approximation. The issue of whether or not the 
runaway mode is significant for snowball growth would 
require an investigation using high resolution coupled N- 
body and hydrodynamic models that exceed the scope of 
this paper. 

In summary, we believe that the interpretation of the 
present-day astero id size distribution as a proof of their 
large initial sizes (jMorbidelli et al.l 120091 ) needs further 
investigation in the light of the results presented here on 
the importance of dust sweeping for early planetesimal 
growth. At the very least, snowball growth opens up the 
possibility that various conditions with different combi- 
nations of R p q and t f could all potentially reproduce the 
asteroids' size-distribution, to the point that the solu- 
tion may be highly de generate. In fact, recent work by 
IWeidenschillind (|2010f ) has already shown that planetes- 
imals with birth-size {R p o) of 0.1 km are also a viable 
initial condition for a model which can well reproduce 
the asteroids' size-distribution. 

Further work would be required, using a detailed nu- 
merical collis ion-and-fragmentationa l m odel (E.g. as im- 
pleme nted in IMorbidelli et al.l (|2009t) or iWeidenschillind 
(|2010T l). to understand whether the observed size distri- 
bution in the asteroid belt can at all constrain the range 
of R p q and tf in the early Solar System. I.e. it may be 
that certain ranges of R p q and t f corresponding to sub- 
regions of the red regime of FigJ5] may give rise to size 
distributions that are better able to reproduce the ob- 
served size distribution in the asteroid belt, while other 
regions can be excluded. 

However, even so, we would still be far from arriving 



at a final answer, as several important questions remain 
unaddressed, including (1) What is the distribution of 
sizes with which planetesimals arc initially created? (2) 
How are the planetesimal births distributed in time and 
space (linear or nonlinear)? (3) To what degree should 
runway growth and gravitational focusing be considered 
in the snowball growth model? All of these factors can 
change the dimension of the red regime of Fig[2] and 
constraining any of these factors will require further de- 
tailed studies of planetesimal formation itself. Last but 
not least, in addition to studies of the formation of aster- 
oids and planetesimals, we may also obtain constraints 
on planetesimal birth sizes through other observations, 
such as those of debris disk aro und A-type and G-type 
stars (|Kenvon fc Bromlevll20Tol ). 

4.1.3. Planetesimal Birth-Rate 

The results shown in Figj2] also allow us to place con- 
straints on the birth rate of planetesimals or their for- 
mation efficiency for each possible formation scenario. In 
addition, they might help overcome some difficulties that 
some of these different formation scenarios encou n ter. 

Recently two new models (jJohansen et al.l 120071 : 
iCuzzi et al.l 120081) have shown that large planetesimals 
can form directly from the concentration of small solid 
pa rticles in the turbulent disk. The numerical simulation 
bv lJohansen et ahl (|2007| ) show that if the feedback of the 
solid particles within the gas is considered then the con- 
centration of solid particles (which is most efficient for 
particles of ~ 50 cm radial size) can grow and be main- 
tained long enough ( the so -called streaming instability, 
lYoudin fc Goodman] (|2005[ )). for the formation of very 
large planetesimals (typically 100-1000 km). Within this 
framework, i.e., with such large initial planetesimals, our 
results (Figf2]) show that the phase of snowball growth 
would be inefficient. In such a scenario, the only way to 
convert a large amount of dust to planetesimals within 1- 
lOMyr is through a very effective planetesimal formation 
process. Nevertheless, one must be careful and note that 
Figf2] ignores the runaway growth mode which might be 
relevant to such large planetesimals. If a runaway mode 
is included, then the snowball timescale, t snow , given in 
Fig(2] should be considered as a upper limit. However, 
as discussed in §4.1.2, the degree to which a runaway 
mode should contribute must rely on future studies with 
high resolution coupled N-body and hydrodynamic mod- 
els. Moreover, one shou ld also note tha t disk turbulence 
is much more intense in Uohansen et al.l (120071) than that 
in Figf2]of this paper. As in Uohansen et al.l (|2007f ) . the 
disk viscosity coefficient is adopted a ~ 10~ 3 , leading to 
H„/H d close to unity according to Eqn.108 of lArmitagd 
(|2010f ) and thus inefficient snowball growth would result. 
It would be interesting to investigate the size and effi- 
ciency with w hich planetesima l s can be formed via the 
mechanism of iJohansen et al.l (|2007|) but with weaker 
turbulence. 

The other large-planetesimal model (jCuzzi "eTaIll2008l) 
also relies on turbul ence which can genera te vortices to 
trap small particles dHeng fc Kenvonll2010D of chondrule- 
size (typically mm-cm). ICuzzi et al.1 ( 20081 ) showed that 
mm-sized particle can be sporadically concentrated by 
disk turbulence, forming large and gravitationally bond 
clumps, which can potentially shrink to solid planetes- 
imals roughly 10-100 km in radius. However, based on 
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the model of ICuzzi et all (120081 ) . iChambersI (j2010l) and 
iCuzzi et al.l (|2010f) investigated the planetesimal forma- 
tion rate in details and found that it should be very low 
in the Solar System: at 1 AU for a typical MMSN disk, 
the time it would take to convert most dust into plan- 
etesimals, tf, would be over ~ 10 10 yr. Snowball growth 
might help in solving this problem by creating a second 
channel by which dust can be converted into planetes- 
imals. Figj2] shows that, for the case of weak turbu- 
lence (this turbulence level in Fig|2] is compatibl e with 
the weak turbulence adopted in ICuzzi et al.l (|2008f )) , the 
time to convert dust into planctcsimals is close to t snow , 
which is much shorter than tf. In fact, t snow can be 
within the disk life-time (1-10 Myr) even for tf > 10 10 
yr. In other words, snowball growth would fully domi- 
nate in terms of transforming dust into planetesimals. 

Beside the two new models described above, there 
are the two traditional models of planetesimal formation 
described in §1. One is the mutual sticking scenari o 
(|Weidenschilling fc Cuzzl 119931 ; iWeidenschillmgl I1997D . 
in which planetesimals form by pair-wise collisions of 
small particles. The other is the instability scenario 
(ISafronovlll9"69l: IGoldreich fc Wardlll973l: lYoudin fcShul 
120021) . which suggests that planetesimals form via grav- 
itational instabilities in a dense particle sub-disk near 
the midplane of the protoplanetary disk. Both scenarios 
remain debated. 

As mentioned in our introduction, the collisional sce- 
nario is challenged by the well-known "meter-barrier" , 
while for the instability scenario it has not yet been es- 
tablished whether a disk of solid particles can become 
dense enough and dynamically cold enough to trigger in- 
stability. While these issues go well beyond the scope 
of this work, the present results show that, for both 
scenarios, snowball growth might represent an interest- 
ing way to allow planetesimal growth to happen despite 
these difficulties. Indeed, in both models planctcsimals 
are thought to form with a radius of 0.1-10 km. For 
this size range, Figj2] shows that snowball growth could 
convert most of the dust mass into planetesimals even if 
formation rates are very low (large tf). This means that 
only a few " lucky seeds" are needed in order for planetes- 
imal growth to proceed. Within the framework of these 
two scenarios, these seeds could either be a few lucky pre- 
planetesimals t hat have overcome t he "meter barrier" for 
some reasons (IBrauer et all 120081: Uohansen et al.l 120081: 
iTeiser fc Wurml 120091: IWettlauferl I2009D or kiiometer- 
sized objects that were formed in some isolated location 
where the critical condition f or gravitational instabilit y 
just happens to be satisfied (jGoldreich fc Ward! I1973D . 
Therefore, snowball growth lowers the requirements for 
those two traditional scenarios: direct planetesimal for- 
mation from dust could be very inefficient and still allow 
planetesimal growth to occur. This hypothesis remains 
to be quantified: future work should focus on deriving the 
planetesimal formation rates within the frame of those 
two traditional scenarios, to see whether they can meet 
these "lowered" requirements. 

4.2. Implication for Planet Formation in Close 
Binaries 

The context in which snowball growth might find its 
most interesting application is that of planet forma- 
tion in binaries. Recent studies have shown that one 



of the major problems for planet formation in close 
binary systems, such as a Centauri, is the intermedi- 
ate stage of planet formation, i.e, the mutual accre- 
tion of km-sized pl anetesimals to form larger planetary 
embryos or cores dThebault et al.1 120061: iHaghighipourl 
120091 : IHaghighipour et al.ll2009|) .The companion's grav- 
itational perturbation, coupled to gas drag, may excite 
large relative velocities between the planetesimals, lead- 
ing to di sruptive collisions wh i ch inhibit their mutual ac- 
cretion dThebault et al.l 120081 120091: IPaardekoooer et al.l 
120081 iMarzari et al.l 120091 : iXie et al.1 l2010T l Recently, 
some mechanisms have been fo und to be somewha t help- 
ful in solving this problem (IXie fc Zhoul [20081 [2001 
IPaardekoooer fc Leinhardt] 120101 ) . but several problems 
remain to be overcome. 

A possible solution to these problems would be 
to let planetesimal collisions only begin when large 
(> 50 — 100 km) ob jects are present in the system. 
iThebault et al.l (|2008[ ) have indeed shown that the grav- 
ity of such large planetesimals is strong enough for them 
to survive the high-speed collisions induced by differen- 
tial gas drag. One obvious way to bypass this problem- 
atic kilometer size range would b e to have planetesimals 
that a re "born big" , be it by the Uohansen et al.l ()2007l ) 
or the[C UZZl "eTalT(l2008l) scenario. However, the issue of 
how these planetesimal formation scenarios might work 
in the very specific context of binaries has not yet been 
addressed. 

Snowball growth offers an attractive alternative way to 
bypass the kilometer size range, regardless of the initial 
size at which planetesimals appear in the disk, because it 
predicts that mutual planetesimal impacts will become 
important only by the time large objects have formed. 
However, it is not possible to directly apply the results 
of the previous sections, derived for single stars, to the 
close binary case. In this section, we address the issue of 
how snowball growth proceeds in the specific context of 
double stars. 

4.2.1. An Example: R p and t h col in a-Centauri 

For the sake of clarity wc consider a disk with a 
lxMMSN surface density and set a = 1 AU, and then 
investigate how snowball growth may occur in binary 
star systems, taking as a representative example the 
case of a Centauri AB, for which the binary separa- 
tion is an = 23.4 AU and the eccentricity es ~ 0.52 
(jPourbaixet al.ll2002l ). 

The main difference between the single star case and 
the binary case is that in the binary case there is a much 
higher relative velocity between the dust and the plan- 
etesimals (v re i). For such a highly eccentric case, we 
can safely neglect the small component of v re i due to 
the differential Kcplerian velocities between the gas and 
the planetesimals (because the gas is pressure supported) 
and basically consider two cases, depending on the un- 
known angle %b between the circumprimary disk plane 
and the binary's orbital plane: 

1. The coplanar binary case, where the binary orbit 
and the gas disk are in the same plane. In this 
case, to a first approximation, the average eccen- 
tricity of planctcsimals (e p ) is equal to their forced 
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eccentricity (5eBa/4as) El, and thus v re i ~ e p Vk ~ 
(5a/4as)e_B«fc, leading to u re ( <~ 840 ms" 1 for the a 
Cen case. This v re i value, and thus also the snow- 
ball growth rate, is roughly one order of magnitude 
higher than in the single star case. We then obtain 
R p ~ 10 -3 km.yr -1 in a weakly turbulent disk with 
H g /H d = 150. 

Note that we here implicitly assume that the gas 
disk is axisymmetric (circular) in the circumpri- 
mary midplanc. This is a rather crude approxima- 
tion as it ignores the reaction of the gas disk to the 
binary pcrturbcr. Nevertheless, as a first order of 
approximation, we believe it is both reasonable and 
convenient for our semi-analytical study in this pa- 
per. In fact, if the reaction of the gas disk is consid- 
ered, v re i would probably be even larger, as shown 
by Paardekooper et al. (2008). Therefore, the v re i, 
as well as the snowball growth rate R p derived in 
our simplified axisymmetric gas disk model should 
be treated as a lower limit. 

2. The inclined binary case, where the binary orbital 
plane is tilted by an angle {is) relative to the 
gas disk plane . In s uch simulated by 

lLarwood et al.l ([1996D . if the Mach number of the 
gas disk is not too high, the gas disk can main- 
tain its structure and undergo a near rigid preces- 
sion. As a planetesimal also undergoes a similar 
precession (regarding the binary orbital plane as 
the reference plane) but with a different rate, the 
relative angle between the gas disk plane and plan- 
etesimal orbital plane will remain approximately 
within the range — 2ig- Taking the medium value 
is as the average, we then have v re i ~ igVk ~ 
520(is/l°) ms . Taking into account this v re i, 
as well as the fraction of time that a planetesimal 
spends moving within the dust-disk, 2Hd/i^aiB, 
then the snowball growth rate for the inclined case 
is R p ~ 10 -5 km.yr -1 , which is independent of both 
is and the turbulence factor H g /Hd- 

Note that, for the inclined case, we have implicitly 
assumed that the vertical excursion of the plan- 
etesimals is higher than the dust disk thickness, 
i.e., i B > 0m(H d /H g )(a/AU) - 25 . However, even 
for a fully turbulent disk (Hd/H g = 1), this condi- 
tion is easily met as long as is > 3°. In addition, 
we also implicitly assume for the inclined case that 
v re i is dominated by is rather than es, which re- 
quires 520(zs/l°) > 840 (i.e. comparing v re i in the 
the inclined and coplanar cases), i.e., is > 1.6° 0. 
Given the above considerations, we shall consider 
that the inclined case discussed here corresponds to 
cases with ig > 3° . Otherwise it should be treated 
as the coplanar case. 

1 In reality, the equilibrium eccentricity can depart from this 
purely dy namical forced value, as it a lso depends on gas drag (see 
Eqn.24 of Paardekooper ct al. (2008)). However, we ignore this 
refinement for the order-of-magnitude feel of the present discussion. 

2 This small value is due to the fact that we are considering a 
region (1AU) close to the primary, where the forced eccentricity, 
which decreases as 1/a, is low whereas the forced inclination, which 
is independent of the semi-major axis, stays at a high value 



The collision timescale, t co i, among planctesimals is 
also different in the binary case. The relation given 
in Eqn llOi which assumes an unperturbed disk with 
equipartition between in and out-of-plane velocities, 
breaks down in binary sys tems. Here we u se the em- 
pirical scaling law given by IXie et al.1 (pOlOl ) (see Eqn.6 
and Eqn.7 in their paper for details) to estimate the plan- 
etesimal collisional timescale in binary star systems as 

t b col ~(0.02 + 2i B )t s col , (14) 

where t s col is the collisional timescale for single star 
systems as given by EqnfTOl 

4.2.2. An Example: The Snowball Phase in a-Centauri 

With the modified R p and t b coll we can derive t snow , 
Rsnow and PMF snow as in the single star case. Results 
are plotted in Fig|3] and Figf?] for the coplanar and in- 
clined binary cases respectively. 

1. Coplanar case. Compared to the single star case of 
Fig|2j we obtain slightly smaller, but still relatively 
large values of PMF snow , i.e., ~ 7% — 50%. R sn0 w 
is little changed, but t snow is about one order of 
magnitude smaller (Fig|3]). These high PMF snow 
and low t snow values imply that snowball growth 
can be efficient, even if the disk life-time in such 
close binaries could be an order o f magnitude 
short e r than in single-s tar systems (jXie fc Zhou! 
120081 : iCieza et al.l 120091 ) . Note that, as for the 
single star case, results depend on the disk tur- 
bulence, which is controlled, to a first approxima- 
tion, by Hg/Hd- Here, we have assumed a value of 
Hg/Hd = 150, corresponding to weak turbulence. 

2. Inclined binary case. PMF snow are here higher 
than in the coplanar case, being close to 50% 
throughout the R p o — tf plane of Fig0] (In fact, 
PMF snow > 50%, but we stop the evolution of the 
system when PMF snow reaches the alarm value 0.5 
as mentioned at the end of section 2.2). Compared 
to the single star case of Fig|2l R snow is smaller 
by ~ 30%, but t snow is larger by about 1 order 
of magnitude (more precisely, by a factor of ~7). 
This implies that the disk lifetime has to be rela- 
tively long in order for snowball growth to be ef- 
ficient. As an example, if planetesimals are born 
small (Rpo < 10 km), then in order for them to 
reach R sno w ~ 50—100 km, i.e., values large enough 
to survive the high-speed collisions induced by the 
companion, then the disk life-time should be no less 
than ~ 10 6 yr. Note that, as R p is independent of 
Hg/Hd (see $4.2.1), the result is here independent 
of the level of disk turbulence. 

The results for both cases considered show that snow- 
ball growth might indeed be an attractive alternative 
to solving the planet- formation-in-binaries dilemma, i.e., 
the high impact velocities that prevent the mutual accre- 
tion of kilometre-sized bodies. It's main appeal is that 
it provides a much safer way for planetesimals to reach 
sizes in the 50-100 km range, after which they are large 
enough to be protected from destructive collisions. The 
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only problem lies with the growth timescales, in partic- 
ular for the inclined binary case for which it might take 
as long as 10 6 yr to reach the safe 100 km size, by which 
time the gas disk might have vanished. However, this 
issue might not be as dramatic as it might appear be- 
cause, even if 50-100 km bodies have not been formed by 
the time the gas disk is dispersed, the enviro nment could 
be fav orable to accretion anyway. Indeed, Kie fc Zhoul 
(|2008|) have shown that, during this gas dispersion phase, 
planetcsimal orbits get re-phased so that, by the time the 
gas has dispersed, relative velocities amongst them are 
very low and accretion friendly. This means that, even 
if snowball growth is switched off (because there is no 
dust sweeping when grains stop to follow gas stream- 
lines) , planetesimal growth can now proceed through the 
classical mutual-accretion channel. 

We note that the 2-D i nvesti gation performed by 
iPaardekooper &; Leinhardtl (|2010T l demonstrated that 
planetesimals can grow to ^50-100 km size with the 
help of dust accretion. Our own investigation extends 
this to the 3-dimensional case, demonstrating that even 
though 3-D collision rates are significantly lower, dust ac- 
cretion (or snowball growth) can still help planetesimals 
to overcome the "km-barrier" for planetesimal accretion 
(or growth) in close binary systems, such as a Centauri. 

4.2.3. limitations 

One possibly problematic assumption of our snow- 
ball growth model regards the efficiency of dust accre- 
tion onto planetesimals, parameterized by Ed defined in 
Eqn0J for which we have implicitly assumed that it is 
not affected by the high v re i values reached in the bi- 
nary e nvironment. Although recent laboratory experi- 
ments (jTeiser fc Wurml I2009H showed that small ( mm- 
sized) particles can be efficiently accreted by large ob- 
jects in a high-speed collision, their test speeds were 
only of a few 10 ms" 1 , which is only relevant for the 
snowball growth in a single star system with v re i ~ 75 
ms -1 . In close binary systems, such as the a Centauri 
shown above, v re i could be as high as 500-5000 ms -1 , 
which is beyond the limit of the test speed in any ex- 
isting laboratory experiment. Therefore, the issue of 
whether dust can be accreted onto small planetesimal 
seeds (i? p o =0.1 — 1 km) for such large relative velocity 
remains unresolved. However, at the very least, it seems 
reasonable to assume that the efficiency of dust accretion 
onto planetesimals will be less affected by high velocities 
than would the efficiency of mutual planetesimal accre- 
tion. 

5. CONCLUSION AND PERSPECTIVES 

We have investigated the transitional phase of planet 
formation, dubbed "snowball growth phase" , in which 
planetesimals are progressively emerging from the dust 
disk but are not yet numerous enough to enter the phase 
of mutual planetcsimal-planctcsimal collision. In this 
snowball phase, the planetesimals move on Keplerian or- 
bits and grow purely through the direct accretion of sub- 
cm sized dust which is entrained with the gas in the pro- 
toplanetary disk and moves with sub-Keplerian velocity. 

Using a simplified model in which the planetesimals are 
progressively created from the dust disk, we find that: 

1. In single-star systems, snowball growth can be a 



significant mechanism if the following conditions 
are met: 

• The dust disk must be thin and dense, i.e., the 
gas disk must be weakly turbulent (see § 13.21) ; 

• The planetesimal creation timcscalc, de- 
fined by the time to convert half the to- 
tal dust mass into planetismals by the 
planctcsimal-formation process alone, must 
be long (roughly tf > 10' yr - see § I3.2[) ; 

• Planetesimals must be born (relatively) small: 
if planetesimals are born larger than 100km, 
then (ignoring the runaway growth mode) 
snowball growth will be inefficient or com- 
pletely absent (see § 13. 3|) . 

// these conditions are met, then we find that 

• Size growth can be significant: 2 orders of 
magnitude in radius (6 orders in mass), before 
mutual planetesimal collisions take over. In 
addition, snowball growth can be fast enough 
for it to happen within the typical (l-10Myr) 
lifetime of a protoplanetary gas disk (see 
§& 

• Snowball growth can be the dominant mode 
for the transformation of dust into planetes- 
imals. Except for highly efficient and fast 
planet-formation scenarios, dust is preferen- 
tially accreted onto already formed kilometre- 
sized bodies rather than consumed for forming 
new "initial" planetesimals (see § 4.1.1). 

• Applying the snowball growth model to a fidu- 
cial MMSN disk for the Solar System and as- 
suming weak turbulence, then the turnover 
from snowball to mutual collisions is likely to 
occur when planetesimals reach a size of 100- 
1000 km, irrespective of their birth size (E.g. 
0.1 - 100km). This result, i.e., that mutual im- 
pacts become significant only when large bod- 
ies have formed, may provide an alternative 
expla nation to the size distr ibution of aster- 
oids (|Morbidelli et all (|2009l) . see §4.1.2). 

• Snowball growth could help overcome some 
still unsolved problems inherent to some 
planet formation scenarios. More precisely, 
it lowers the requirements on the planetesi- 
mal formation mechanism itself, by allowing 
planetesimals to grow, and dust to be trans- 
ferred onto them, even if the mechanism by 
which they are formed is very inefficient (see 
§4.1.3). 

2. For close binaries, snowball growth offers a way to 
overcome one major problem for planet formation 
in binaries, i.e., the high impact velocities that pre- 
vent kilometer-sized planetesimals from accreting 
each other. Indeed, this dynamically excited envi- 
ronment, which is hostile to mutual planetesimal 
accretion, is on the contrary favorable to growth 
by dust sweeping. If efficient enough, snowball 
growth allows planetesimals to grow large enough, 
50-100 km, to be protected from mutually destruc- 
tive impacts. Two main cases can be distinguished: 
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• In the coplanar binary case (is < 3°), snow- 
ball growth is at its most efficient (a little 
more efficient than the single star case, if one 
assumes the same level of disk turbulence) but 
remains sensitive to disk turbulence, i.e. weak 
turbulence is still required for efficient snow- 
ball growth (see FigEJand §4.2.2). 

• In the case of an inclined binary (ig > 3°), 
snowball growth can still be efficient, although 
slightly less efficient than in the case of a sin- 
gle star with weak disk turbulence(i? 3 /i/d = 
150), with the added advantage that it is in- 
dependent of disk turbulence (see FigU] and 
§4.2.2). 

Let us stress again that our model is highly simplified 
and that the present work, whose main objective was 
to identify a new and potentially significant mechanism, 
should be considered as a first step towards more de- 
tailed studies. In future work, the modeling of snowball 
growth will need to be improved in a number of aspects. 
The most important improvement will be to have a self 
consistent model of the planetcsimal size evolution under 
the coupled effect of snowball growth and mutual accre- 
tion, instead of the simple analytical growth law assumed 
here. This requires the development of a global statis- 



tical particle-in - a-box m odel, in the spirit of t hose of 
IWeidenschilimgl (fl997l ) or iMorbidelli et all (|2009l) . span- 
ning a wide range in sizes and incorporating all pos- 
sible outcomes for dust-planetesimal and planetesimal- 
planetesimal collisions. Other improvements which could 
be made to the model include, (i) the role of gravitational 
focusing of dust onto planetesimals, in particular to what 
extent it is (or is not) hampered by the action of gas, and 
(ii) the changes to Ed arising from the physics of high ve- 
locity impacts of dust onto large bodies, especially in the 
specific case of close binaries where v re i may reach values 
of a few 100m. s" 1 . 
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Weak turbulence (H g /H d =150) Intermediate turbulence (H g /H d =15) Strong turbulence (H g /H d =1.5) 
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Fig. 1. — Planetesimal mass fraction (PMF), mass- weighted average radii {(Hp)) and planctesimal collisional timescale (t co i) as a function 
of time for the cases with weak (left), intermediate (middle) and strong (right) turbulence. Two subcases with tf = 10 5 yr (in green, 
denoted with a triangle) and tf = 10 11 yr (in red, denoted with a square) are plotted in each panel. The dash-dot lines denote the location 
of time = i co ; line, so that the vertical coordinates of the triangles (or squares) indicate the location of t Bnow , Rsnow and PM F snow . N.B. 
the red and green lines in the 3 horizontal-middle panels which should be absolutely overlapped are deliberately plotted with a small offset. 
We see that increased tf (i.e., low planetesimal formation efficiency) allows more snowball growth to occur before collisions commence, 
increasing R snovJ by ~ 1 order of magnitude. In addition, we note the effect of strong turbulence in decreasing the efficiency of snowball 
growth, resulting in lower P M F snow and Rsnow occurring before planctcsimal-planctcsimal collisions commence. 
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Fig. 2 . — t snow (solid contours filled with colors) and R S now (dashed contours) in the RpQ — tf plane at 1 AU for 1 MMSN with 
Hg/Hd = 150 (weak turbulence) in a single star system. The i? p o-coordinate is schematically divided into three size ranges for the initial 
planetesimals, forming via traditional models (Goldreich & Ward 1973; Weidenschilling &; Cuzzi 1993) and via two new models which 
involve turbulence (Johanscn ct al. 2007; Cuzzi et al. 2008). Snowball growth is most efficient in the top-left corner of the plot, where 
planetesimals of initial size 0.1 km can grow via direct dust accretion (snowball growth) to ~ 100 km within 10 6 yr, before planetcsimal- 
planetcsimal collisions commence. Note that, for t s now in the ~ 10 6 - 10 7 yr range (the typical life-time of a protoplanetary disk), we 
obtain 100 < Rsnow < 1000 km regardless of the initial value R p o- 
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Fig. 3. — Same as Figf2]but for the coplanar binary case (is = 0) and Hg/Hj = 150. The binary system parameters are those of a 
Centauri. Rsnow is close to its value for the single star case FigEJ but t snow is about one order of magnitude smaller, implying that 
snowball growth is efficient even if the disk life-time in close binary systems is an order of magnitude shorter than in single-star systems. 
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Fig. 4. — Same as Fig[3]but for the inclined binary case with is = 10°. The binary system parameters are those of a Centauri. Note 
that the present results are independent of the turbulence factor Hg/H^. Compared to the single star case of Fig[2] Rsnow is smaller by 
~ 30%, but t sn ow is larger by about 1 order of magnitude. If, for example, one expects planetesimals to be born small (R p q < 10 km) but 
want them to have a large Rsnow (10-100 km) in order that they can survive high-speed collisions, then the disk life-time should be no less 
than ~ 10 6 yr. 



APPENDIX 

MASS-WEIGHTED AVERAGE RADII: (R P ) 

As we assume a linear planctcsimal formation rate (see EqnQ} and a linear growth rate (see Eqn[7|). thus, at a given 
time t, the distribution of the planetesimals' radii will be given by R p o, RpO + dr, R p q + 2 x dr, R p q + i x dr, R p o + 

(n — 1) x dr, where dr is the radii increase during a time interval of tf.indi i- e -j dr — R v tj.ind, and n denotes the total 
number of planetesimals. The total mass of these planetesimals is 

M tot = -jP* zJ^ + * x dr ) 3 

= -^-P*(R 3 p0 n + 3R 2 p0 dr ^ i + 3R pQ dr 2 V i 2 + dr 3 V i 3 ). (Al) 
We define the mass- weighted average radii, (R p ), which satisfies 

M tot = —p*{R p ) 3 n. (A2) 
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TABLE 1 
List of Notions. 



Variables 


Meaning 


Definition 


O.B 


Scmimajor axis of the orbit of binary system 


First paragraph of §4.2.1 




Ice boundary (a semimajor axis of ice condensation) 


Just behind Eqn.fSJ 


/'•-; 


Dust accretion efficiency factor 


Eqn.lflJ 


e-B 


Eccentricity of the orbit of binary system 


First paragraph of § 4.2.1 




Planctesimal average orbital inclinatio 


Second paragraph of §4.2.1 


fd 


Enhancement factor of from a nominal disk (MMSN) 


Just behind Eqn.fSJ 


fice 


Enhancement factor of due to ice condensation 


Just behind Eqn.fSJ 


//,. 


Scale height of the dust-disk 


Eqn.fTJ 


//, 


Scale height of the gas-disk 


Eqn.fSJ 


i-B 


Inclination of the orbit of binary system 


Second paragraph of §4.2.1 




Planctesimal average orbital inclinatio 


Eqn.ffljJ 


M p0 


Initial individual mass of the planctcsimals 


Eqn.lfTJ 


M p 


Planctesimal mass growth rate via dust accretion 


Eqn.fjJ 


M p 


Planctesimal mass 


Eqn.fjJ 


N 


Surface density of the planctcsimals 


Eqn.lfTJ 


n 


Planctesimal volume density 


Eqn.fTuTl 


PMF 


Planctesimal mass fraction in the total solid mass 


Eqn.lTJJ 


PMF snow 


PMF at the end of snowball phase 


Eqn.lfTJ]) 


Rpo 


Initial individual radial size of the planctcsimals 


Third paragraph of § 2.1 


Rp 


Planctesimal radii 


Eqn.lflJ 


Rp 


Planctesimal radial growth rate via dust accretion 


Eqn.fSJ 


{R P ) 


Mass- weighted average radii of the already- formed planctcsimals 


Eqn.f8j 




{Rp) at the end of snowball phase 


Eqn.fJJJ 




Timescalc for the onset of planctcsimal-planctcsimal collision 


Eqn.fTuTl 


t b , 

COL 


t co i in binary star systems 


Eqn.fTUl 


L col 


t co i in single star systems 


Eqn.fTUl 




Formation timescale of individual planctcsimals 


Second paragraph of §2.1 


'/ 


Characteristic planctesimal formation timescalc 


Eqn.lfJJ 




Duration of snowball phase 


Eon.fTTll 




Transition timescale from dust-dominated to planctesimal dominated 


First paragraph of § 4.1.1 




Planctesimal escape velocity 


Second paragraph of §4.1.2 


Vk 


Local Kcplcrian velocity 


Eqn.lO 


Vrel 


Relative velocity between planctesimal and dust 


Eqn.fTJ 


AV 


Planctesimal average relative velocity 


Eqn.fTuTl 


£ p 


Planctesimal formation efficiency factor 


Eqn.lfTJ 


P* 


Planctesimal internal density 


Eqn.fSJ 


s d 


Dust surface density 


Eqn.lfTJ 



Note. — The table list is on the alphabetical order of the initial of the variable name. Several variables in Greek are listed at the end 
of the table. 



Comparing Eqn IAll and Eqn lA2[ we have 

= iKo + ?^o(« - I)* 1 + R P o(n ~ l)(n - \)dr 2 + -n{n - lfdr 3 ] 1 / 3 . (A3) 



l p0 ~ r, 1 l p0 

Considering 



l)dr ~ (n — -^)dr ~ ndr ~ R p t, 



then EqnlA3l can be rewritten as 



3 -1 

(R p ) ~ [i? 3 + -Rlo(Rpt) + R p o(R p t) 2 + - (Rpt) 3 ] 1 / 3 
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R 



po 



1/3 



Rpt 



(A4) 



